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We report on the preparation and detection of entangled states between an electron spin 1 /2 and 
a nuclear spin 1/2 in a molecular single crystal. These were created by applying pulses at ESR (9.5 
GHz) and NMR (28 MHz) frequencies. Entanglement was detected by using a special entanglement 
detector sequence based on a unitary back transformation including phase rotation. 
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The entanglement between two spins 1 /2 is at the heart 
of quantum mechanics. Ever since a so-called "paradox" 
was formulated by Einstein, Podolsky and Rosen (EPR) 
101 , referring to local measurements performed on the in- 
dividual spins of a delocalized entangled pair, proper- 
ties of entanglement and its consequences for quantum 
physics has been discussed in great detail [||. In the 
context of quantum information processing (QIP) entan- 
glement has been considered as a resource for quantum 
parallelism (speedup of quantum computing) || and 
quantum cryptography 0. A number of these quan- 
tum algorithms have been demonstrated in NMR (nu- 
clear magnetic resonance) quantum computing ^, ^ . 

In this contribution we report on the experimental 
preparation and observation of the entangled states of 
an electron spin S = 1/2 and a nuclear spin / = 1/2 
in a crystalline solid. The spins considered here are 
a proton and a radical (unpaired electron spin) pro- 
duced by x-ray irradiation of a malonic acid single crys- 
tal ||ll| . This leads to the partial conversion of the CII2 
group of the malonic acid molecule to the radical *CH 
where the dot marks the electron spin. In a strong 
magnetic field the following four Zeeman product states 
\msmj) = I TT), I Ti), I it), I ii) exist where the 
arrows label the ±1/2 states of the electron and the nu- 
clear spin. Equivalently we will use a qubit labelling as 
\msmi) = |00), |01), |10), |11). The energy level dia- 
gram corresponding to the electron-proton spin system is 
shown in fig. |l|, where we have also indicated the possible 
ESR {Ams = ±1) and NMR transitions (Am/ = ±1) of 
the individual spins by solid arrows. 

What we are aiming at are states of the type 
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These represent all four possible entangled states of a two 
qubit system, also called the Bell states of two spins 1/2. 
They correspond to a superposition of the states in fig. |l| 
connected by dashed arrows. 

Electron spin resonance (ESR) was performed at X- 
band (9.49 GHz) at T = 40 K. The low temperature 
was chosen only for reasons of signal-to-noise ratio. The 
two well resolved ESR lines due to the 'CH proton de- 
pend on the orientation of the single crystal and were 
observed at magnetic fields of 338.2 mT and 339.2 mT 
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FIG. 1: Schematic diagram of the four energy levels of a two 
spin system with 5 = 1/2 and / = 1/2. The solid arrows 
denote allowed transitions. The dotted arrows indicate for- 
bidden transitions, corresponding to entangled states. The 
phase dependence of the quantum states under z-rotations is 
also indicated (see text). 



with a linewidth of 0.5 mT for the ESR and about 1 MHz 
for the ENDOR (electron nuclear double resonance) line. 
This orientation corresponds to a principal axis of the 
hyperfine tensor. There are two different proton NMR 
transitions. We applied pulsed ENDOR techniques to 
one of them at the frequency of 28.05 MHz. 

In the high temperature approximation we express the 
Boltzmann spin density matrix as /5b = (1 — KB)ji + 
Kb ■ Pp with Kb = /iB^o/^B?^ (for g — 2) and where 
the pseudo Boltzmann density matrix is defined as pp = 
{jl — ■^Sz) which corresponds to an equal polulation of 
the states | |t) and | H) {p^ = P4 = 1/2) and equiva- 
lently I Ti) and I TT) {Pi = P2 = 0) with tr{pp} = 1. 
We used here the assumption that the Larmor frequency 
of the nuclear spin loqi is much smaller than the Lar- 
mor frequency of the electron spin ujqs- Note that 
the pseudo-pure density matrix poo can be expressed as 
Poo = ii + + \iz + Sjz = |00) (00| corresponding to 
the pure state ^ = |00). 

In what follows we will prepare density matrices corre- 
sponding to the Bell states according to eqn. (|l|). Since 
we will apply selective transitions we need to consider 
only a three level subsystem for each of the Bell states. In 
order to prepare all four states and $^ we only need 
to apply transition selective excitations to either of the 
three level subsystems 1,2,3; 1,2,4; 1,3,4 or 2,3,4. Here 
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we restrict ourselves to the sublevels 1,2,4 for creating 
the <I>^ states and 1,2,3 to create the vE'''' states. As an 
example we treat the state in detail. It corresponds 
to the well known Einstein Podolsky Rosen (EPR) state 
The preparation of '5^ proceeds by first applying a 
selective 7r-pulse to the 1 <-+ 3 transition to create the 
following pseudo pure populations: pi = 1/2, p2 = 0, 
P3 = of the 1,2,3 three level subsystem. The corre- 
sponding pseudo pure density matrix of the three-level 
subsystem (1,2,3) represents the pseudo pure state |tT)- 
The creation of the states can now be achieved by the 
pulse sequence P/^(=F7r/2) followed by Pg'^{—TT) where 
the plus sign in P/^(=F7r/2) creates the VP" state. Here 
we use the abbreviation j{f3) for a selective pulse at 
the transition j ^ k with rotation angle (3. The label 5* 
refers to an electrons spin transition, whereas the label / 
stands for a nuclear spin transition. The corresponding 
unitary transformation results in 

I tt) "^'^-^'^ ^(1 TT> + 1 Ti» "^'^"^ ^(1 Ti) - I iT». 

(2) 

In order to create '^'^ the plus sign in P/'^(±7r) must 
be chosen. For completeness we note that the states 
can be created by following the same line of reasoning 
when starting from the sublevels 1,2,4 with a prepara- 
tory TT-pulse at the 2^4 transition. Next the pulse 
sequence Pj^(±7r/2) followed by P|'*(— tt) is applied cre- 
ating the <I>^ state except an overall minus sign. Other 
scenarios using the other sublevels are possible and will 
be presented in a more extensive publication. 

In order to prove that the state has indeed been 
created we apply a density matrix tomography which is 
based on the phase dependence of the entangled state 
as already sketched in fig. |^. The phase factors noted 
there represent the phase dependence of the correspond- 
ing states under rotation about the quantization axis (z- 
axis). This corresponds to the unitary transformations 
Us^ — exp(— 1015^) and C//^ = exp(— i(/)2/z) leading to 
Us^UiJmsmi) = exp(-i(0ims + (j)2mi))\msmi) . A 
single ESR transition (Ams = ±1) will have a phase 
dependence 4>i under z-axis rotation, whereas a single 
NMR transition (Am/ = ±1) will have a phase depen- 
dence 02- Each of the entangled states and is 
characterized by a linear combination 01 ± 02 of both 
phases. This is another manifestation of the fact, that 
these states are global states and no local measurement 
on the single qubits reveals any information about the 
entangled state. 

Since the entangled state is not directly observable we 
need to transform it to an observable state. Our en- 
tangled state detector therefore corresponds to a uni- 
tary back transformation comprised of e.g. Pg^(— vr) fol- 
lowed by P/^(— 7r/2) applied to the state 4"". In order 
to distinguish entangled states from other superposition 
states we encode their characters in the phase depen- 
dence under z-rotation as discussed before. Therefore we 



apply phase shifted pulses of the type Pg^{—'K, 0i) and 
Pj^(— 7r/2, 02) which corresponds to a rotation about 
the z-axis by the angles 0i and 02. The complete ^'^ 
detector sequence now reads 

;7*(0i, 02) =Pf (-V2, 02)Pi'(-7r, 0i). (3) 

Since our observable is the ESR transition intensity, de- 
tected via an electron spin echo, the entangled state to- 
mography corresponds to the evaluation of the following 
signal strength for 'I'^ 

5*(0i, 02) = tr {2Sf (If IvE--) {^-\Up] (4) 

where S],^ is the fictitious spin 1/2 of the 1^3 transi- 
tion. With the current definitions we obtain the expres- 
sion 

5d*('/'l,02)-^(l-COs(01-02)). (5) 

Alternatively we have used for the detection of the <I>^ 
states the sequence 

J7*(0i, 02) = Pf(-V2, 02)P|"(-7r, 0i) (6) 

leading to a detector signal for $+ 

^*(0i, 02) = tr {2^24 ^1 i^+^^^+i fj^^Y 

with 

^l(</'l,</'2) = ^(l-COs(0i-f02)). (8) 

A more detailed discussion of the phase dependence of 
the detector signal will be presented in a more extended 
publication. 

The phase shifts were implemented by incrementing 
the phase of the individual detection pulses in consecu- 
tive experiments according to 0j (n) = Acjj n At with j = 
1,2. The artificial phase frequencies Awj = 27rAi/j were 
arbitrary chosen as i^i = 2.0 MHz and V2 = 1-5 MHz. Ex- 
amples of different phase increments are shown in fig. |^. 

Four different sets of phase variations were chosen to 
demonstrate the individual and combined phase frequen- 
cies. In fig. H a, b we have set (a) vi — Q oi (b) V2 = Q 'vci 
order to demonstrate the 02 and 0i dependence as a ref- 
erence. The corresponding spectra (see fig. ||a, b) are ob- 
tained after Fourier transformation. These would also be 
observed for non-entangled superposition states of either 
ESR [vi) or NMR transitions {v2)- The characteristics 
of the entangled states shows up in the combined phase 
dependence (see eqns. and This is demonstrated 
for the VE"" state in fi-g.^ c where the interferogram al- 
ready shows the phase difference behavior which is even 
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FIG. 2: Phase interferograms versus time n At for four dif- 
ferent sets of experiments, (a) <j)i = 0, (b) (^2 = 0, (c) 
(1)1 ^ 0, (t)2 ^ Q for the entangled state ^~ with phase depen- 
dence 101 — 02| (see eqn. (d) 01 7^ 0, 02 7^ for entangled 
state with phase dependence |0i -|- 02 1 (see eqn. |^ 
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FIG. 3: Fourier transform of the phase interferograms shown 
in fig. ^ for phase frequencies vi = 2.0 MHz and U2 = 1.5 MHz 
(see text), (a) 1^2 = 1.5 MHz, (b) i/i = 2.0 MHz, (c) ji/i -1/2I, 

(d) \U1 +U2\. 

more clearly evident from the spectral representation of 
fig. I^c showing a line at the difference frequency vi — V2- 
In a similar way the states were created which under 
the corresponding tomography sequence show a phase 
dependence as 0i + 02 which is evident from the interfer- 
ogram (fig. |^) and more clearly from the phase spectrum 
(fig. leading to a spectral line at vi + V2- 



Here we have used an ensemble of electron- nuclear spin 
pairs with a mixed state density matrix. Like in all so 
far published NMR quantum computing experiments the 
corresponding entangled state would be better termed 
pseudo entangled states |^, However, we point 

out that the same pulse sequences could be applied to 
a pure state electron-nuclear spin pair in order to create 
the discussed entangled states. Also the same Bell state 
detection sequences proposed here would apply. In fact 
the same phase dependence would be observed for pure 
states. Moreover, the experiments presented here would 
reach the quantum limit (see Warren et al. |]l2j) ac- 
cording to the PPT (positive partial transpose) criterion 
]l3| , at tanh(/ki;5/2fcBT) = \/2 which corresponds to 
a temperature Tq = 2.576 K for an ESR frequency of 
95 GHz. 

Theoretically the scenario reported here might seem to 
be obvious, however, experimentally this type of experi- 
ment is rather demanding because of the extreme differ- 
ence of the ESR and NMR linewidth and the duration 
of the microwave and radio-frequency pulses. We have 
used the following typical values: 32 ns for the ESR and 
1.6 /j,s for the NMR 7r-pulses. Entanglement is achieved 
with these in 832 ns. However, because of the broad ESR 
and NMR lines a complete excitation of the transitions 
could not be obtained with these pulses. This leads to er- 
rors in the effective rotation angles of the corresponding 
pulses. As a consequence of this an incomplete creation 
of the entangled states results, leading to residual sepa- 
rable states. These are expected to vary as 0i and 02- 
In order to estimate this effect we calculate the resulting 
phase dependence of the detector signal (electron spin 
echo) for the deviation 5j of the 7r/2-pulses at ESR {5i) 
and NMR {52) frequencies which results in 

'S'd(0i, </'2) = ao + ai cos 01 -I- 02 cos 02 

-t-ai2cos(0i - 02) (9) 

with oi = — i(Si(l — 5^)52^ 02 = —j5i52 and 012 — 
- 62) up to fourth order. 

We note that the pulse errors introduce cos0j (j=l,2) 
dependencies in addition to the expected cos(0i — 02) 
dependence of the entangled state. In order to reduce 
the cos (jjj dependencies in the detector signal we have 
applied a phase cycling procedure where the two signals 
Sd{<t>ii 4>2) and 5^(01 -I- TT, 02+7r) were added. This elimi- 
nates the cos dependencies but retains the cos(0i — 02) 
dependence of the entangled state. This procedure was 
applied to the data presented in figs. ^ and |^. 
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